Dirac Point Structure in a Bose-Einstein Condensate in a Honeycomb Optical Lattice 
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We study the Bose-Einstein condensate in a honeycomb optical lattice within Bogoliubov theory 
and find that for a k = condensate, the Dirac points appear in the Bogoliubov excitation spectrum 
when < /3 < 2, which illustrates that the bose-bose interaction does not change the Dirac point 
structure but only give a modification of the velocity of the Dirac cone. When the bosons are 
driven to condense at k = K, however, we find that the topology of the Dirac points will be altered 
by arbitrary weak interaction. Furthermore, we find that the next-nearest-neighbor hopping in an 
isotropic and an anisotropic lattice has different effects to the dynamics of the condensate and it 
should be taken into account when the lattice is not sufficiently deep. 
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In recent years, the Dirac point structure has arisen a 
lot of interests in condensed matter physics[ - t], non- 
linear optics [4-6] and cold atomic physics [7, 8]. The 
Dirac point structure corresponds to many important 
phenomena in physics, such as the room-temperature 
quantum Hall effect in graphene[ ], topological edge state 
in topological insulator [ ], and conical diffraction in 
honeycomb photonic lattices [ ]. To realize such a struc- 
ture in a cold atomic system, several groups proposed 
to load fermions on different optical lattices [ - ] and 
recently, in experiment, L. Tarruell et al has realized it 
in a tunable honeycomb lattice [15]. Before this realiza- 
tion of Dirac point structure with Fermi gas, the hon- 
eycomb lattice was first investigated with Bose-Einstein 
condensates [16, 17], although new quantum phases were 
observed, no signatures of Dirac points were observed. 

In this paper, we study the Bose condensates in a hon- 
eycomb optical lattice and find that in the tight-banding 
limit, for a k = condensate, stable Dirac points appear 
in the two lowest bands even in the presence of the inter- 
action if the anisotropy j3 (the meaning of j3 is given in 
the following) is within the region (0, 2). However, if the 
bosons are condensed at k = K, where K is the momen- 
tum where the two bands touch, even within Bogoliubov 
theory, we find that no matter how weak the interaction 
is, the Dirac points will be altered by the interaction, 
which agrees with the conclusions obtained in Refs.[6, 18] 
by numerical methods. Furthermore, we find that for 
a lattice that is isotropic but not sufficiently deep, the 
next-nearest-neighbor (NNN) hopping effectively affects 
the dynamics of the condensate but nearly does not af- 
fect the dynamics around the Dirac points. However, 
for a anisotropic lattice, the NNN hopping will alter the 
topology of the Dirac points and make the Dirac points 
no longer well defined. Therefore, to observe the Dirac 
points, it's better to use an isotropic lattice. 

Model-TYie honeycomb optical lattice, which consists 
of two sublattices A and B, can be realized by three de- 



tuned standing-wave lasers, with the optical potential 
given by[ ] 

V{x,y)= ^ FiSin^ [/cl(x cos l9i sin l9i) + 7r/2] , 

i=l,2,3 

where Oi = 7r/3, O2 = 27r/3, Os = 0, and /c^ is the optical 
wave vector. With different Vi, the honeycomb lattice 
can be either isotropic or anisotropic[ ]. In this work, we 
consider single-component bosonic atoms in this lattice. 
For bosons, the intra-species collisions is dominated by 
s-wave scattering. In the following, we consider tq << 
k~^as^ where ro is the effective interaction length, k = 
^ is the wave vector, and is the s-wave scattering 
length. Furthermore, we first assume the lattice is deep, 
therefore, the effective Hamiltonian is given as the Bose- 
Hubbard model 
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where < ij > represents the nearest neighbor (NN) sites, 
di and bi denote the bosonic mode operators for the sub- 
lattices A and B, respectively, /i is the chemical potential 
and U describe the on-site interaction between bosons. 
The tunnelling rates tij, in general, depend on the tun- 
nelling directions in an anisotropic honeycomb lattice. 
In this paper, both the isotropic case and the anisotropic 
case are analyzed in detail. 

First, we make a Fourier transformation for the Hamil- 
tonian. With dl = {l/VN)^^exp{ik ' Ai), = 
{l/^/N)^^exp{i]^ • Bi)bj^, where N is the number of 
sites of the sublattice A (or B). For a strong optical lat- 
tice, the restriction to the lowest Bloch band is justified. 



and the Hamiltonian can be given as 



(2) 



where = Ep ^p-q^P = «^^)- ^{^) = 

ELi^^exp(ik-bs), with bi = (l/V3,l)(a/2), b2 = 
(1/^3, -l)(a/2), and bs = (-a/V3,0), a is the lattice 
spacing. In this paper, we set ti = t2 = t and = j3t. 
Furthermore, we set t and a as the energy unit and length 
unit, respectively. 

For the free case, i.e. /7 = 0, this Hamiltonian can be 
rewritten in a 2 x 2 matrix form. 



-M 0(k) 

/.*(k) -^Ji 



by making a transformation ((ik,/5k)''^ 
where the matrix V takes the form 



V 



1 



i'^(k)i -m 



T^(ak,&k)^, 



(3) 



the Hamiltonian is diagonalized and the energy eigenval- 
ues are 



e±{k) = -^i±\m\- 



(4) 



Since the bosons are condensed into the zero crystal 
momentum, S-(0) should take zero, and therefore ji = 
— 1 0(0)1 . Such a non-zero value of chemical potential is 
common in lattice models and it guarantees the energy 
of excitations to be positive, which is necessary for the 
condensates to be stable. 

The dispersion relations are determined by |(/>(k)|, 
which is directly related to the hopping amplitude and 
takes the form 



2cos(fc,)+4/3cos(^)cos(^). 




(5) 



From Eqs.(4), (5) and Fig.l, we see, for < /3 < 2, the 
two bands touch at several points and Dirac points ap- 
pears. In particular for the isotropic case where j3 =, the 
dispersion has the same Dirac structure as the graphene 
material (as shown in Fig.l). The Dirac points are purely 
due to the structure of the honeycomb lattice, without 
any relation to the quantum statistics. For /3 > 2, the 
two bands no longer touch with each other and they are 
gapped. 

When the interaction is tuned on, the chemical poten- 
tial has a shift: fi = — |0(O)| —\(j){0) \ -\-Un^ where n is 
the density of the system. In the Bogoliubov approxima- 
tion and neglecting the constant condensation energy, the 



Hamiltonian can be written as Hbg = iXlk^k^k^k, 
where 5^ = (aj^, a_k, ^k, ^-k) and 



/eo^U U (/)(k) \ 

U so^U 0(k) 

0*(k) so^U U 

\ (/)*(k) U so^U J 



Here we assume the filling factor of bosons is unit and 
So = |(/>(0)| just for convenience. Based on this Hamilto- 
nian, we find that the dispersion relations are 



£±(k) = yjel + 2Ueo + ^(k)|2 ± 2{eo + C/)|^(k)|. (6) 

If we take U to be zero in Eq.(6), it is same as Eq.(4). 
Furthermore, for < /3 < 2, the two bands keep touching 
and the Dirac point structure is not destroyed by the 
interaction, which is shown explicitly in the following. 

For < /3 < 2, there are two regions which we are 
interested in. The first one is around the k = 0, by 
expanding the momentum k around (0, 0) and up to the 
second order of kx and ky, the dispersion relation (6) 
becomes 



Vyky, 



(7) 



where = ^JWp/2{2^ p) and Vy = ^JuJ2. The sound 
velocity is anisotropic and dependent on the anisotropy 
that is characterized by j3 of the lattice. The second 
region is around the Dirac points. Following the same 
procedure, we obtain 



^±(q) 
^±(q) 



So ± 



vlql 



- VyQy, 



(8) 



■ 2Uso ± 



Sq + U 
V4 + 2C^£o 



(9) 



where q is a small momentum away from the Dirac 
points K = {kl.k^), i.e. {k^.ky) = [k^ + q^.k^ + qy). 

Vx = and Vy = ^^T^^^^^T^- The anisotropy of 

the velocity is a consequence of the anisotropy of the lat- 
tice. Removing the constant energy, we see that s±{q) 




FIG. 1: (Color online) (Right) the isotropic case, i.e. (5=1. 
The two bands touch at six points, with two different kind 
of Dirac points. (Left) the anisotropic case with /3 = 3, the 
bands are gapped. 
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represents the standard energy-momentum relation for 
the relativistic Dirac particles, with Vx and Vy replacing 
the light velocity. Therefore, there is real Dirac struc- 
ture around the touching points. Furthermore, compared 
Eq.(9) to Eq.(8), we see that the interaction has only 
modified the velocity of quasiparticles around the Dirac 
points. 

Before discussing how to detect the Dirac point struc- 
ture in experiments, we have to realize that here the 
physics is in the ultracold and dilute region, the quan- 
tum depletion is very small, e.g. for ^^Na at a typical 
density of 10^^ cm^, the quantum depletion is 0.2%. The 
quantum depletion for two dimensional system is even 
more smaller [21]. Therefore, unlike other cases, here the 
physics is dominated by the condensates and low energy 
excitation. As a result, the effects of the Dirac points 
are very weak for a k = condensate to be detected. 
This is the reason why no signatures of Dirac points is 
observed in experiments [16, 17]. To observe the effects 
of the Dirac points, we have to drive the bosons to con- 
dense at k = K. However, when the bosons are driven to 
condense at k = K, there is modulational instability, this 
can be confirmed by the discrete nonlinear Schrodinger 
equation (DNLS)[^:2], 



ih 



dt 

dipB 



^-gf- = ^ts^pAS + Ul'tpBfll'B, 



(10) 



where s{s) denotes the nearest-neighbor vectors from 
A(B) to B(A). The stationary solutions of Eq.(lO) are 



plane waves ipQ^'^expli^k • A(B) — i^t)], of frequency v 



±|0(k)|+[/ (here we assume ng = 



A\2 



B\2 



no = 1). To check the stability of these states, we perturb 
the carrier wave with small amplitude phonons: = 
(V^^'^+7/A,Be^^-A(^)+^i,^e-^^-A(^))e^(k-A(^)--*),The 
DNLS excitation spectrum is determined by 



\ T,q B^HqEq, where = (aj,, a_q, 6^, 6_q) and 
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The energy spectrum is given as 

e^{q) = V(|^(q + K)|±2C/)|0(q + K)|. (12) 

We see no matter how weak the interaction is, imaginary 
always appears in the long- wavelength limit, i.e. q ^ 
0. As a result, the topology of the Dirac points will be 
altered by the long- wavelength excitations, which agrees 
with the result obtained in Refs.[6, 18]. All these results 
indicate that to observe the Dirac point structure in Bose 
condensates is a challenge. Recently, a way to handle this 
problem of stability was introduced [ ] and the vortex 
patterns suggested there can be used as a probe of the 
Dirac point structure in experiments. 

The study above is based on the tight-binding model. 
We know the tight-binding model is based on the as- 
sumption that the lattice is deep enough to guarantee 
that the Wannier wave function is well-localized. How- 
ever, for a lattice that is not sufficiently deep, the tight- 
binding model with just the NN hopping is much less ac- 
curate than the one with both NN hopping and the NNN 
hopping[ ], and therefore, the NNN hopping should be 
taken into account for such not sufficiently deep lattice. 
In the following, we will consider the NNN hopping to 
see what effects the NNN hopping has to the condensate. 

After we introduce the NNN hopping into Eq.(l), the 
excitation spectrums take the form 



e±(k) = -/2-Mk)±|0(k)|, 



(13) 



- [2^^ + mi. + I0(k + q)P + |^(k - q)p]A2 
-2iU + Ck)[|0(k + q)P - |<^(k - q)|2]A + i^l + 2/7^)' 
-U'm + q)* + </.(k - q)|2 + |0(k + q)|2|0(k - q)|2 
-(e^ + 2C/a)[|^(k + q)|2 + |0(k-q)|2] 



0, 



(11) 



where (fk = |(/>(k)|. When k = 0, the spectrum is given as 
A± = ^yel + 2Ueo + |(/>(k)|2 ± 2(^0 + /7)|(/)(k)|, the same 
as Eq.(6) and therefore it is just the Bogoliubov spec- 
trum. When k = K, imaginary will naturally appear in 
the excitation spectrum, and as a result, the condensate 
is modulationally unstable. 

In fact, when the bosons are condensed at k = K, we 
can even see that the Dirac points are destroyed by the 
interaction within the Bogoliubov theory. Based on the 
Bogohubov theory, the Hamiltonian is given as Hbg = 








'-•4 


2 y4\6 



FIG. 2: (Color online) (a) 1,2,3 stand for three different kind 
of hopping, (b) The common parameters for both the solid 
line (black) and the dashed line (red) are /3 = X = 1, U = 0.5 
and /cx = 0, while 7 = for the solid line and 7 = 0.1 for the 
dashed line. 
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where 

ji = -27(1 + 2A) -/3-2, 
/i(k) = 27 (cos ky + 2 A cos cos 

where 7 is the NNN hopping constant in unit of t and 
A is the anisotropic of the NNN hopping. Based on the 
anisotropy of the lattice and the anisotropy of the NN 
hopping, the hopping amphtudes along the path labelled 
by 2 and 3 in Fig. 2 (a) are equal but different to the one 
along the path labelled by 1. From Eq.(13), it is direct 
to see that the two bands still touch. However, when we 
expand the dispersion relation around the Dirac points 
and keep only the terms related to the lowest order of 
momentum, we find 

£±(q) = -27(/3 - X)sjl -^qy ± ^^^+^y, (M) 

this indicates that if the anisotropy of the NNN hopping 
is different from the NN hopping, the Dirac point struc- 
ture is no longer well defined even though the two bands 
still touch with each other (as shown in Fig. 3). In the 
harmonic approximation, we calculate the NN hopping 
and the NNN hopping and find when /3 7^ 1, A is not 
equal to j3 (if we further consider the next-next-nearest- 
neighbor (NNNN) hopping, the NNNN hopping will in- 
duce a small gap for the case of {3 ^ 1). Therefore, to 
observe the Dirac point structure, it's better to choose 
the system to be isotropic, i.e. j3 = \ = 1. In the follow- 
ing, we set /3 = A = 1. 

When the interaction is tuned on, the dispersion rela- 
tions change into 

£±(k) = ^el + 2/7ek + |(/>(k)|2 ± 2(ek + t/)|(/>(k)|, (15) 

where ek = — /i — /i(k). By expanding the dispersion 
relation around /c = 0, we obtain 



e_(k) = ^vlkl + vlkl, (16) 




FIG. 3: (Color online) (Right) The isotropic case, i.e. (5 — 
A = 1; (Left) The anisotropic case, i.e. (3 ^ \. 



we see the sound velocity is effectively increased due to 
the NNN hopping (as illustrated in Fig. 2(b)). By ex- 
panding the dispersion relation around the Dirac points, 
we obtain, 

.±(q) = A(7, U) ± -^^^^2^^T^, (17) 

where ^(7, U) = 9(1 + 7)^ + 6/7(1 + 7) and ^(7, U) = 
3(1 + 7) + /7. By a direct calculation with concrete pa- 
rameters, we find that unlike the case around k = 0, the 
NNN hopping nearly does not affect the velocity around 
the Dirac points (also shown in Fig. 2(b)). Based on these 
results, we obtain the conclusion that for an isotropic lat- 
tice, the NNN hopping mainly affects the dynamics of the 
modes around the k = 0, it nearly does not affect the dy- 
namics of the modes around the Dirac points; However, 
for an anisotropic lattice, the NNN hopping not only af- 
fects the dynamics of the modes around the k = 0, but 
also alters the topology of the Dirac points and makes 
the Dirac points no longer well defined. 

From the above analysis, we see that to guarantee that 
the tight-binding model gives an accurate description to 
the condensates, the lattice indeed has to be sufficiently 
deep to guarantee that the NNN hopping is far smaller 
than the NN hopping, i.e. 7 <C 1. 

In summary, we have studied the Bose-Einstein con- 
densates in a honeycomb lattice and found that for a 
k = condensate, the Dirac points appear in the Bo- 
goliubov excitation spectrum when < /3 < 2, and when 
the bosons are condensed at k = K, the topology of the 
Dirac points will be altered by arbitrary weak interac- 
tion. Furthermore, we find that the NNN hopping of 
not small strength has obvious effects to the dynamics 
of the condensates, and therefore it should be taken into 
account when the lattice is not sufficiently deep. 
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